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UNIQUENESS ISSUES FOR EVOLUTION EQUATIONS WITH DENSITY 

CONSTRAINTS 

SIMONE DI MARINO AND ALPAR RICHArD MESZAROS 


Abstract. In this paper we present some basic uniqueness results for evolutive equations under 
density constraints. First, we develop a rigorous proof of a well-known result (among specialists) 
in the case where the spontaneous velocity field satisfies a monotonicity assumption: we prove 
the uniqueness of a solution for first order systems modeling crowd motion with hard congestion 
effects, introduced recently by Maury et al. The monotonicity of the velocity field implies that 
the 2—Wasserstein distance along two solutions is A-contractive, which in particular implies 
uniqueness. In the case of diffusive models, we prove the uniqueness of a solution passing 
through the dual equation, where we use some well-known parabolic estimates to conclude 
an contraction property. In this case, by the regularization effect of the non-degenerate 
diffusion, the result follows even if the given velocity field is only I/°° as in the standard Fokker- 
Planck equation. 


1. Introduction and preliminaries 

Recently, modeling crowd behavior has received a lot of attention in applied mathematics. 
These models actually are in the heart of many other ones coming from biology (for instance cell 
migration, tumor growth, pattern formations in animal populations, etc.), particle physics and 
economics (see for example the recently introduced models of Mean Field Games, [21, 22, 23]). 
For more details on these models we direct the reader to the non-exhaustive list of works [5, 6, 
7, 8, 11, 12, 13, 14, 15]. In all these models the question of congestion can play a crucial role. 
Indeed, from the modeling point of view one could have some ‘singularities’ if individuals want 
to occupy the same spot. In this paper, we will work with equations which model some type of 
congestion effects in crowd motion models (for a more detailed description of these models we 
direct the reader to the references [25, 27, 31]). These systems read as 

dtpt - + V • (ptVj) = 0 

< Pt < 1, p \ t=Q = Po 
A’adm (pt) ■ 

In the above system pt represents the density of a crowd (at time t) that moves in C for 
a time T > 0 accordingly to the prescribed velocity Ut, a given field that everyone would follow 
in the absence of the others, u € {0,1} is just a parameter: if u = 0, one deals with first order 
systems (the density of the population is just transported by some vector field) introduced in 
[25], while for u = 1 one has a second order system (in addition to the transport, the population 
is also affected by some randomness, modeled by a non-degenerate diffusion) studied in [31]. In 
order to preserve a density constraint (we suppose that pt does not exceed a given threshold, 
let us say 1), at each moment one modifies to vj, a field that is the closest to Ut (in the 
sense) and it is an admissible velocity field, i.e. the set adm(pf) represents the fields that do not 
increase the density on already saturated zones {pt = 1}. 
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In generic congested models a pressure is also acting (Darcy’s law): this pressure is preventing 
congestion and, according to various models, it is an increasing function of the density. In the 
usual soft congestion models we have some power law of the type p = p™ (porous medium 
equation). However this choice cannot prevent p to be very high, while it is clear that a crowd 
of people cannot have a density higher than a certain threshold (studies say that the maximum 
density is 4.5 people/m^), which we put for convenience equal to 1: we will refer to the constraint 
p < 1 as an hard congestion effect. We shall present later how the pressure field appears in the 
models considered in this paper. 

Introducing general multivalued monotone operators /3, the above systems can be written in 
a compact form as 


dtPt + V • (utpt) = Apt in (0, T) x H 

(1.1) < {utpt — Vp) • h = 0 on (0,T) X dQ 

Pt G /3(pt). 

The two cases i/ = 0 and u = \ correspond to two special operators /3i and (52, namely 


h{p) 


/ 

0 

< [0,+oo] 
+ 00 

\ 


if p < 1 
if p = 1 
if p > 1, 




/ 

P 

< [l,+oo] 
+ 00 

\ 


if p < 1 
if p = 1 
if p > 1; 


notice in particular that /32(p) = /3i(p) + p. The hard congestion effect is due to the fact 
that (5i{p) = Too whenever p > 1: in fact, this will force p to be always less than 1. It 
is worthwhile noticing that this problem has some features in common with the Hele-Shaw 
model and the Stefan problem (see for instance [17, 18]), namely in both problems there is a 
degenerate monotone operator linking the density and the diffusive part. However there are also 
big differences: in our case we treat also the convection term, which can depend not smoothly on 
the position, while usually in the Hele-Shaw models there is a source term for the mass, which 
we do not treat here since we want to model a crowd moving inside a domain H. This modeling 
assumption is also the reason why we consider the no-flux Neumann boundary condition, since 
we want neither people exiting nor entering the domain. 

Moreover our equations can be seen also as a quasilinear elliptic-parabolic system with very 
degenerate nonlinearity (see [33, 16, 4]), for which the issue of uniqueness is still an open problem 
when we add a driving vector held without imposing an entropy condition. In this context we 
may write our equations using the variable u = p + p as 


dtg{u) + V • (4>(t, X, u)) = Abi{u), 

where bi{u) = {u — 1)+, g{u) = u — bi{u), b 2 {u) = u and ^{t,x,u) = \it{x)g{u). In the case 
i = 1 we have a double degeneracy since both bi and g have a hat part, instead for the case 
i = 2 we have only a degenerate part in g{u). However again the (possibly) rough coefficient 
in $ rules out results already present in the literature. For problems of this form we expect an 
contraction result, see for example [4]; notice also that we do not need a concept of entropic 
solution since the equation reduces to a linear one in the joint variable (p,p) (see the systems 

(1.2) below). Thus the usual concept of weak solutions can be considered for our purposes. 

In fact, we will be dealing always with the following reformulations of the systems, using the 
fact that p G (5{p) if and only if p < 1 and p G 77^(12) (see its dehnition in (1.3)): 

(1.2) 

dtPt + V • (utpt) = Apt in (0, T) x Q 
< {utpt — Vp) ■ h = 0 on (0, T) X (9H < 

^Pt < 1, Pt G 


dtPt + V • (utpt) = A{pt + Pt) 
(utpt - Vp - Vpt) • h = 0 
Pt < 1, Pt G 


in (0, T) X Q 
on (0, T) X dQ 
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In Subsection 1.1 we will derive in another way these systems: that derivation justifies also 
the regularity assumption of pt- 

A very powerful tool to attack these type of macroscopic hard-congestion problems - where we 
impose a density constraint on the density of the population - is the theory of optimal transport 
(see [38, 37]), as we can see in the recent works [25, 26, 27, 31, Ij. In this framework, the density 
of the agents satisfies a continuity, or a Fokker-Planck equation (with a velocity field taking into 
account the congestion effects) and can be seen as a curve in the Wasserstein space. 

Our aim in this paper is to prove some basic results of uniqueness of the solutions in this 
setting. As far as we are aware of, this question is a missing puzzle in full generality in the 
models studied in [25, 26, 27, 31]. Let us remark that the uniqueness question is crucial if one 
wants to include this type of models into a larger system and one aims to show existence results 
by fixed point methods, as it is done for Mean Field Games in general for instance (for example 
as in [35]). 

We will treat two different cases, with very different approaches: in the first one we simply 
consider a crowd driven by a given velocity field and subject to a density constraint. In this case 
the assumption that the velocity field is monotone will be crucial in order to prove a A-contraction 
result for the solutions, that will imply uniqueness. In the second case we add a diffusive term, 
which models some randomness in the crowd movement (see [31] for recent developments and 
existence results in this setting); in this case we prove a standard contraction property 
passing to the dual problem and proving there existence for sufficiently generic data. In this 
case a major role is played by the regularizing effect of the Laplacian, which allows us to prove 
the result even if the velocity field is merely bounded. 

Let us underline that the core of both methods is classical in the literature. 

(i) On one hand we use the differentiation of the squared Wasserstein distance along two 

solutions of continuity-type equations, and then use a Gronwall-type argument to show 
a contractive property and hence uniqueness of solutions for these evolution equations. 
This type of proof is very common in the optimal transport theory; nevertheless our 
analysis requires some finer nontrivial new ideas, mainly because of the appearance of 
the new pressure variable which we can formally identify as a term associated to the 
subdifferential of the indicatrix function of the (geodesically convex) set {p € : 

p < I}. Moreover, we believe that by the nature of our problem the machinery of the 

contraction through the doubling and re-doubling of the variables (successfully used 
for instance in [33] and [4]) seems to be very difficult, complicated and heavy to adapt 
to our situation, while the method proposed by us is simple and elegant. In addition, in 
this context we give new proofs for some of the used results from the theory of optimal 
transport. 

(ii) On the other hand, in the second order case (when we add a non-degenerate diffusion 
term into the model) the idea to pass to the dual equation to show the uniqueness is very 
similar to the techniques developed already in [9] for instance; the method to obtain the 

contraction follows the same lines of [3]. However the PDF studied here seems to 
be new in this context and the result of uniqueness is per se interesting in the theory of 
crowd motion. 

We underline moreover that the two methods used in the two models are mutually exclusive: 
the W 2 distance along two solutions of the second order model would be contractive if the vector 
field would be monotone, but we do not require this assumption; the parabolic estimates used for 
the second model highly rely on the fact that one has a non-degenerate diffusion in the system, 
which is clearly not the case for the first order model. 

We remark also that we expect a contraction result also in the first case, since it can be 
seen as a doubly degenerate quasilinear elliptic-parabolic equation [33]. In Section 4 we provide 
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a sketch of this result in the case in which the velocity field is monotone: we underline that we 
use the uniqueness proved in Section 2. 


1.1. The density constraint: admissible velocities and pressures. In order to model 
crowd movement in the macroscopic setting with hard congestion, we work in a convex bounded 
domain C such that |11| > 1. The evolution of the crowd will be analyzed by the evolution 
of its density, which is assumed to be a probability measure on Q. The condition we want to 
impose is a bound on the density of the crowd, which is considered to be always less than or 
equal to 1. In particular the set of admissible measures will be denoted by 


/Cl := {p G : p < 1 a.e.} 


(here and after we identify p with its density with respect to the Lebesgue measure 

As for the velocities, we need to impose that the density is not increasing when it is saturated: 
informally we would say that v is an admissible velocity for the measure p G /Ci if V • v > 0 in 
the set {p = 1} and v • h < 0 on dVl where h is the outward normal. In order to make a rigorous 
dehnition we have to introduce the set of pressures: 

(1.3) Hl{n) := [p G H^{n) : p > 0, p(l - p) = 0 a.e.}. 

Then, using the integration by parts formula, formally for p G Hp{Q,) and v admissible one 
should have 

0 < / pV • vdx — / pv • = — / V • Vpdx 

Jq JdU Ju 

and so we can define 


adm(p) := < V G 



• Vp dx < 0 for every p G 



In the sequel we will denote by T’adm(p) the L^(II,M'^)—projection onto the cone adm(p). 

Now, in order to preserve the constraint p < 1, we impose that the velocity always belongs to 
adm(p) and so a generic evolution equation with density constraint will be 


dtPt + V • (pivt) = 0 
Pt < 1, vt G adm(pt). 


One of the simplest such model is when we have a prescribed time-dependent velocity field 
and we want to impose that the velocity is the nearest possible to u*, time by time. This 
describes a situation where the crowd wants to have the velocity but it cannot, because of 
the density constraint, and so it adapts its velocity trying to deviate as little as possible: this 
will result in an highly nonlocal and discontinuous effect. The first order problem hence reads 
as 


(1.4) 


dtPt + V • (ptVi) = 0 in (0, T) x fl 
< Pt < 1, p|t=o = Po 

-Padm(pt) [Ut] j 


where the first equation is meant in the weak sense and the minimal hypothesis in order to 
have a well dehned projection is u G L^([0,T] x II). In the following lemma we characterize the 
projection of the velocity field: 


Lemma 1.1. Let p G 1^(11) such that p < 1 a.e. and let u G L^(II;M'^). Then there exists 
p G Hp{Q) such that T’adm(p) ['*^] = u — Vp. Furthermore p is characterized by 
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Proof. Let us set K = {Vp : p G Hp{Q)}. It is easy to see that K is a closed cone in 
Let us recall that the polar cone to K is dehned as 

;= |v G L2(0;M‘^) : J v ■ Vqdx < 0, y Vq £ K 

By the definition of the admissible velocities we have adm(/9) = K°. Moreau decomposition 
applied to K and K° gives 


u = Pk[u\ + Pko[u] Vu G 

This proves the claims (i) and (ii), since these are precisely the conditions of being the projection 
onto a cone. □ 


Corollary 1.2. Lemma 1.1 (i) implies that 
in particular we get 



|Vppdx+ / \u—Vp\'^dx, and so 


[ |Vpp dx < f 

Ju Jq 


u\^ dx, 


/ 


-Padm(p) [’*^11 dx < 


/ 


up dx. 


Using Lemma 1.1 we get that if (p, v) is a solution to (1.5) then — Vp*. Now, using 

that pt'^pt = Vpi, we have that (p,p) is a solution to 


dtPt + V • (ptut) = Apt, in (0,T) x U 
Pt < 1, Pt G 


Thus we found again equation (1.2), where we imposed also the no-flux boundary condition, 
which is the one expected when we model a crowd in a closed room Ll. Now we are ready to 
give the formal definition of solution to our problem: 

Definition 1.1. Let u G L^([0, T] x Pi) and let po G /Ci. Then we define a solution to (1.5) to 
he a couple (p,p) G L°°([0,r] x O) n ^C([0, T]; ( '^(O), IT's)) x L‘^{[d,T];H^iPl)) with p\t=o = pq 
and such that: 

• for all (j) G C')N(M'^) and for all 0 < r < s <T we have 


n {ut{x) — Vp{x)) ■ \/(j){x)pt{x) dx dt = / (l){x)ps{x) dx — / cj){x)pr{x) dx] 

\ J Q, J Q, 


• we have 0 < p < 1 /or j2f^L[0,T]^Af'^Lll-a.e. (t,x) andpt G Hp^{Q) /or Af^L[0,r]-a.e. 

t. 


Remark 1.2. We remark that the boundary condition implied by our definition is homogeneous 
Neumann. Note that we assume for the density p to be also an absolutely continuous curve in the 
space of probability measures equipped with the Wasserstein distance^. Like this it is meaningful 
to assume p\t=o = Po- 

In [27, 25] and [36] the following regularity hypotheses have been assumed to show the exis¬ 
tence result: u G or n = —VD for a A— convex potential D and in both cases no dependence 
on time. 


^see the subsection on optimal transport below and [2] for further details. 
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1.2. The diffusive case. Recently (see [31]) a second order model for crowd motion has been 
proposed. In a nutshell, it consists of adding a non-degenerate diffusion to the movement and 
imposing the density constraint. This leads to a modified Fokker-Planck equation and, with the 
notations previously introduced, it reads as 


( 1 . 6 ) 


dtpt - ^Pt + V • (ptvt) = 0 

< /Ot < 1, p\t=Q = Pq 

= -fadm(pt) [Ut], 


where is - as before - the desired given velocity field of the crowd. Introducing the pressure 
gradient in the characterization of the projection, one can write system ( 1 . 6 ) as 


dtPt + V • (ptut) = A(pt pt) 
Pt < 1, p\t=o = Po, Pt G 


where, as before we equip the equation with the natural no-flux boundary condition on cffl; the 
rigorous definition of solution is similar to that of Definition 1.1. We notice also the fact that 
we have in particular pt + pt^ h{pt) and so we are in fact solving equation ( 1 . 1 ) with /3 = (32■ 
Under the assumption that D is convex and u E L°°([0,T] x D) it has been shown (see [31, 
Theorem 3.1]) that the system (1.7) admits a solution {p,p) £ T] x D) x T^([0, T];H^{^1)). 

In addition [0, T] ^ t pt is an absolutely continuous curve in the 2-Wasserstein space (see 
Subsection 1.3). We direct the reader to [31] for further details on this model. 


1.3. Optimal trasport. Here we collect some facts about optimal transport that will be needed 
in the sequel. Given Xi and X2 two measurable spaces and T a measurable map between them, 
we say that the Borel measure p2 is the push forward of the Borel measure pi through T and 
we write p2 = T^Pi, if P 2 {A) = piiT~^{A)) for every measurable set A C X 2 . 

Given two measures p € i^{Xi) and u G ^(X 2 ), we define n(p, u) as the set of 7 € i^{Xi x 
X 2 ) such that (vri)(j 7 = p and ( 712)07 = where TTj is the projection to the i-th coordinate: 
these measures are called transport plans between p and u. 

A particular example of transport plans is given by the transport maps: whenever we have 
a measurable T : Xi ^ X2 such that T^p = 12 we have that the induced plan jt = (id, T)^p 
belongs to n(/i, zz). Let us summarize some well-known facts about optimal transport in the 
following theorem (see for instance [38] or [37]). 

Theorem 1.3. Let D C 6 e an open bounded set and let p,it £ 3 P{Q). Let us consider the 
following quantities 

(P) A(/i, i/) = inf I / \x-y\^d-f : -f £U{p,u) 


(D) 

B{p,u) = sup\ / (p{x)dp+ / tf{x)di' : £ Cb{Ll),(p{x) + ftiy) <-\x - y\‘^ 

IJn Jn ^ 

where C'b(D) denotes the space of bounded continuous functions on Q. We will call (P) the 
primal and (D) the dual problem. 

(i) There exists at least a minimizer for the primal problem (the set of minimizers is denoted 

by Ilo{p,r')) and there exists also a maximizer (y?, V’) dual problem. 

(ii) A{p,u) = 2B{p,i') and we will call W^ip^r') the common value. 

(iii) We can choose a maximizer {ip, tp) of (D) such that p and are Lipschitz in Ll and also 

such that ^\x\^ — p{x) is a convex function on the convex hull ofQ. If p then its 

gradient is a map T{x) = x — Xp{x) such that T^p = v and whose associated plan is the 
unique optimal plan, that is Ilo{p,v) = {yr}- 
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(iv) If^l is convex, ^(0) endowed with the Wasserstein distance W 2 is a geodesic space and 
if fi the geodesic between g, and v is unique and it is described by 

[0, l]3t<-^g.t-= (id + t{T - id))^fj,. 

The space of probability measures equipped with the Wasserstein distance W 2 is called the 
Wasserstein space. We denote it by W 2 := (^(n), W 2 ). 


Structurally, the following two sections contain our main results. Section 2 is devoted to the 
uniqueness question for first order models where the main tool is the theory of optimal transport. 
In Section 3 we investigate the uniqueness issue for second order models, using PDE techniques. 
Finally, in Section 4 we discuss an approach that could lead to an contraction result in the 
first order case as well. 


2. Monotone vector fields in the first order case 

Let C be a bounded convex domain. In this section we suppose that the desired velocity 
field u : [0, T] X n — )■ of the crowd is a monotone vector field in L^([0, T] x II; M'^), i.e. the 
following assumption is fulfilled: there exists A € M such that for all t € [0, T] there exists a 
Borel measurable set At ^ Q (possibly depending on t) with full measure, i.e. \ At) = 0 

and 


(HI) {\it{x) - \it{y)) ■ {x - y) < X\x - y\^, x,y e At. 

The following contractivity results are not very surprising in the Wasserstein context. In 
practice we show that (Lemma 2.1) 

{-Vp : 

where the indicatrix function of /Ci and 9^2 denotes the Wasserstein subdifferential (see 
[2]): then we exploit the geodesic convexity of in order to get the contraction properties. 
However, in order to let the reader understand clearly the proofs we will omit the W 2 technical 
language. It would be interesting to adapt these tools to more general Hele-Shaw problems. 

Although a first written version of these results is essentially contained in [30] (Section 4.3.1), 
here we simplified and clarified some of the proofs. A key observation is the following lemma 
(see also Lemma 4.3.13 in [30]): 

Lemma 2.1. Let Ll be a convex bounded domain 0 /and let po,Pi G ^{Ll) two absolutely 
continuous measures such that po < 1 pi < 1 a.e. Take a Kantorovich potential p from po 
to Pi and p G Then 

/ V(p-Vpdx= / V(p'Vpdpo>0 . 

J Q Jo. 

To prove this result we consider the following extra lemma: 

Lemma 2.2. Let LI be a convex bounded domain of and let po,pi G .^(H) two absolutely 
continuous measures such that po < 1 Pi < 1 o..e. Take a Kantorovich potential <p from po 
to Pi and p G H^{Ll). Let [0, 1] 3 1pt be the geodesic connecting po to pi, with respect to the 
2-Wasserstein distance W 2 . Then we have that 


_d 

dt 




Vp • Vpdpo. 
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Proof. We know (using the interpolation introduced by R. McCann, see [29] or Theorem 1.3 
(iv)) that pt = {x — t'Vip{x))^po for all t G [0,1] and so we have 


_d 

dt 


{*= 0 } 


/ 

Jn 


pdpt = lim / 


p{x — tS/(p{x)) — p{x) 


dpQ{x) 


= — lim 
t —^0 


n 
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Vp{x — sS/ip{x)) ■ Vip{x) ds dpo{x) 


= - lim / At{Vp) ■ VipdpQ{x), 

Jn 

where the second equality is easy to prove, for fixed t, by approximation via smooth functions 
and, for t G [0,1], we denoted by At : the linear operator 

rt 


1 r 

At{h){x) = - / h(x — sV<y9(a:)) ds. 
^ Jo 


Now as a general fact we will prove that At{h) h strongly in as t —>■ 0, for every 

h G First of all it is easy to see that ||^t|| < 1: indeed 

1 


/ 

Jn 


\At{h)\^dpo< 



In Jo 

fl 


\h{x — sV<y9(x))f dsd/3o(a;) 


d/?s(x)ds< / \h\‘^dx. 


Here we used the fact that since po,pi < 1 a.e we have also pt < 1 a.e. for all t G [0,1]. 
Now it is sufficient to note that for every e > 0 there exists a Lipschitz function such that 
\\he-hh2 < s, and so we have 

\\At{h) - h\\L2^ < \\At{h — h^)\\L],^ + ||h - < 2e + tL||V(p||j;^2^, 

where L is the Lipschitz constant of h^. Taking now the limit as t goes to 0 we obtain 
hmsupj_,.o \\At{h) — h\\i 2 ^ < 2e; by the arbitrariness of e > 0 we conclude. 

Now it is easy to finish the proof, since Vp G and so 


_d 

dt 


{i=0} ' 


pdpt = - lim / Ht(Vp) • Vpdpo = - / Vp-Vpdpo- 
! Jn 


/ 

Jn 


□ 


Proof of Lemma 2.1. Let [0,1] 3 t i-A pt be the Wasserstein geodesic between po and pi. We 
know that pt < 1 a.e. for all t G [0,1] and in particular it is true that 

/ pdpt < / pdx = / pdpo, 

J Q, Jq J q. 

which means that the function [0,1] 3 t i—>■ / pdpt has a local maximum in t = 0, hence its 

Jn 

derivative at 0 is non-positive. 


Given this, the claim is a consequence of Lemma 2.2. 


□ 


We will also need a regularity lemma on the continuity equation: by definition a curve pt 
satisfies a continuity equation with velocity if for every cp G C'q(M'^)^ the application [0, T] 3 


^we recall that for us C'o(R'^) is the closure of C^(R‘*) with respect to the norm 
equivalently, the set of 0 £ C^(R'^) such that lim|j.|_>oo l<?i(®)l + |N(/)(a;)| —>■ 0 


+ ||Vfli||oo or, 
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t (->■ / dpt is absolutely continuous and its derivative is / v* • V(/>dpi for almost every time 

jRd. ijjd 

t. We will prove that under some integrability assumption there exists a universal full-measure 
set of differentiability, even when (j) £ 

Lemma 2.3. Let pt he a weakly continuous curve of probability measures on let us sup¬ 
pose that the continuity equation dtpt + ^ • {ptVt) = 0 holds with a velocity field Vt such that 

/ / dptdt < -boo. Then there exists a set T C [0,T] such that J^^([0,T] \ T) = 0 and 


/ 


( 2 . 1 ) 


lim — 
h 


(fdpt+h- / (pdpt]= S/(f)-Vtdpt Vt 

jR'i / JR'i 


£ T 


for every £ C'q(M'^). Moreover if we further have that pt < 1 a.e. and / / \vt\'^ dpt dt < -boo 

Jo Jr^- 

then we can require that (2.1) also holds for every (f> £ H^{ 


Proof. Let us prove the following general statement: for a given separable Banach space B and 

a curve x* £ L^{[0,T]; B*), there exists T C [0, T] such that T is a set of Lebesgue points for 

the map t e->- x^{b) for every b € B, and moreover .if^([0,T] \ T) = 0. This can be proven 

easily by choosing a dense subset {bn) C B and then taking as the set of Lebesgue points of 

t !->• x^{bn) and Tq as the Lebesgue points of t i—)• Hx^H, and then we take T = 

every 6 £ B and e > 0 let us consider i £ N such that ||6j — 6|| < e and then taking t £ T we 

have 

^ Jt 5 ^ ^ lksl|dsj-b^ ^ \xl{bi) - x*t{bi)\ds. 

Now, taking the limit as <5 —>■ 0 and using the properties of T we get 

limsup— / \x*{b)-x*{b)\ds <2e\\x*\\ 

5^0 2.0 Jt_s 


and by the arbitrariness of e we conclude that t is a Lebesgue point for x*{h). Now it is easy to 
conclude thanks to the fact that we know that 


f 4 >dpt+h- [ (/)dpt= f f Vfi-Vsdps, 

JRd jRd Jt jRd 

and noticing that x* ■. j ■ Vg dps satisfies the assumption x* £ L^([0, T]-,B*), when B 

jRd 


is the Banach space C'q(M'^) in the first case and H^{M.'^) in the second case. Notice that if we 
follow the construction of T for the case B = CQ(ffi‘^), this set of times will work also for H^{W^) 
since Cq is dense in H^. □ 


Now we are in position to prove the main theorem of this section, namely: 


Theorem 2.4. Suppose Q C is a bounded convex domain, u is a vector field satisfying 
Assumption (HI) and let po £ /Ci be an admissible initial density. Let us suppose that there 
exist (p^,p^), (p^,p^) two solutions to the system 


( 2 . 2 ) 


dtpt + V • {pt{^k - Vpt)) = 0 in (0,T) x Ll 
< Pt < 1, Pf > 0, (1 - pt)pt = 0 a.e. in [0,T] x H, 
, ^ 14=0 = Po, 


p* £ L‘2{[0,T]; H^{Ll)) and pi £ i^{Q) for i £ {1,2}, where the first equation is supposed to be 
satisfied in duality with C^{W^) (see Definition 1.1) in order to take into account the boundary 
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conditions. Then, and a.e. In particular, under the same assumptions, we can 

say that there exists a unique pair {p,v) that solves (1.4). 


Proof. We associate to the two curves p] and pf a continuity equation (1.5) with the correspond¬ 


ing vector fields vj = u^—Vp} and = uf — Vpf. Let us compute and estimate -^-W^ipjjpf); 

we refer also to [2, Theorem 8.4.7] for a more general statement, but we prefer to include a sim¬ 
pler proof in this case. We know that t i-A W^ipjjPt) is absolutely continuous: let us consider 
a time t for which its derivative exists and also such that t G n T^, where T* is the set for 
which (2.1) is satisfied for the continuity equation for p*. Then we know that, for all s G [0,r] 
we have 


(2.3) 


[ Ptdpl+ ! iptdpl, 
^ Jn Jn 


where {ipti'f’t) is a pair of Kantorovich potentials for pi and p^. In particular we have equality 
in (2.3) for t = s and so, since both sides are differentiable by hypothesis, their derivatives are 
equal. Hence we get 


where we used the fact that, since Q is bounded, we can assume (ft and ift to be Lipschitz and 
in particular they belong to which allows to use Lemma 2.3. 

We also know that thanks to Theorem 1.3 there is a pair of optimal transport maps Tt{x) = 
X - V(pt(x) and St{y) = y - Viptiu) such that {Tt)^pl = pi and [Stl^pl = p] for all t G [0,r]; 
moreover Tt is the inverse of St (in the appropriate almost everywhere sense). Using this, we 
have VijjtiTtix)) = —V(pt(x) and in particular using the change of variable formula y = Tt{x) 
we get 

[ Viftiv) • Ui(j/) dpl{y) = - [ S/ipt{x)- ut{Tt{x)) dpl{x). 

Jn J S7 

We can use this to split the formula for the derivative of W^ipl, pl)/2. We use the result of 
Lemma 2.1 and then rewrite the term regarding in terms of transport maps and see: 

= f vl-Viptdpl+ f vf-S/ilJtdpl 

= [ Uf Vpt dp] + [ Uf V'lft dpi - [ Vpj ■ Vpt dp] - [ Vpl ■ Vift dpi 
J Q, J Q. J Q. 

< f Vpt(x) • [ut(x) - Ut(ri(x))] dpi 

Jn 

< f {x-Tt{x))- [ui(a:) - ut{Tt(x))] dpi 

Jn 

<\f \x- Ttix)\'^dpl < XWi{pl,pl). 

Jn 

Gronwall’s lemma implies that 


WUpIpD < c^^^WUpIpI). 

Since p^ = Pq = po a.e., the above property implies that p^ = p^ a.e. in [0, T] x 12. From this 
fact we can easily deduce that A{pl — pi) = 0, for a.e. t G [0, T] in the sense of distributions. In 
particular pi — pi is analytic in the interior of Q. Moreover, both pi and pi vanish a.e. in the 
set {pI < 1} which has a positive Lebesgue measure greater than |n| — 1 > 0. Thus, p^ = p^ a.e. 
in [0, T] X n. The claim follows. □ 
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Remark 2.1. The existence result for system (2.2) was obtained in different settings in the 
literature. On the one hand, if u = —VD (for a reasonably regular potential D), the existence and 
uniqueness of a pair {p,p) can be obtained by gradient flow techniques in W2{0) (see [2, 25, 36] 
under the assumption that u is monotone that translates into a convexity for D. On the other 
hand, if u is a general field with regularity, the existence result is proven with the help of a 
well-chosen splitting algorithm (see [27, 36] j. 

Nevertheless, combining the techniques developed in [31] on the one hand, and the well-known 
DiPerna-Lions-Ambrosio theory on the other hand, we expect to obtain existence result for (2.2) 
for more general vector fields with merely Sobolev regularity and suitable divergence bounds. 

Remark 2.2. The monotonicity assumption (HI) is not surprising in this setting. We remark 
that the same assumption was required in [32] to prove contraction properties for a general class 
of transport costs along the solution of the Fokker-Planck equation in 

dtp-■ (Bp) = 0, p\t=o = po, 

where the velocity field R : —>■ was supposed to satisfy the monotonicity property (HI). 

We note also the fact that we can allow for moving domains Qt (considering always a no-flux 
boundary condition): in fact our proof never uses that the domains are fixed but uses only 
convexity at any fixed time. This generalization has been used in [10] for proving uniqueness for 
an evolution equation with density constraint driven only by the boundary of the moving sets. 


3. Bounded vector field in the diffusive case 


We use Hilbert space techniques (similarly to the one developed in [9, 3, 34]; see also Section 
3.1. from [35]) to study the uniqueness of a solution of the diffusive crowd motion model with 
density constraints described in Subsection 1.2 (see also [31]). Moreover we can expect that 
this holds under more general assumptions in the presence of a non-degenerate diffusion in the 
model. 

Let u : [0, T] X ^ be a given vector held, which represents again the desired velocity 
held of the crowd, 11 C a bounded open set with boundary, po € .^(H) the initial density 
of the population such that 0 < po < 1 a-®- iu ^ aud let us consider the following problem 

r - AptV • (Padm(pt)[ut]pt) = 0, in(0,T)xll; 

(3.1) 

L plt=0 = Po, 0 < Pt < 1, a.e. m H, 


equipped with the natural no hux boundary condition. Introducing the pressure variable, equiv¬ 
alently the above system can be written as 


dtpt - Apt - Apt + V • (utpt) = 0, in (0, T) x H, 

(3.2) p|t=o = Po, in H 

(Vpt + Vpt — Utpt) • h = 0, on dTt, for a.e. t € [0,T], 
for a pressure held pt G It has been shown in [31] that under the assumption that 

(H2) u G L°°([0,r] X 11;M'^) 


the systems (3.1) and (3.2) have a solution. More precisely there exist an absolutely continuous 
curve [0,T] 9 t i-A pt G W2 and pt G for a.e. t G [0, T] (in particular p G L°°{[0,T] x H) 

and p G L^([0, T]; R^(ll)) such that (p, p) solves (3.2) in weak sense (see (3.3)). 

Our aim in this section is to show the following theorem: 
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Theorem 3.1. Let u satisfy (H2). Then there exists a unique pair {p,p) G L°°{[0,T] x LI) x 
L‘^{[0,T]] (Ll)) that solves (3.2) in the weak sense (3.3). Moreover for every solution {p^,p^) 

and {p^-,p^) to (3.2) we have the — contraction property 

[ \p^(.T,x) - p‘^{T,x)\dx < [ \p^{d,x) - p^{d,x)\dx. 

Jo. J Q. 

Proof. The existence of a solution {p,p) was obtained in [31]. Now let us show the uniqueness of 
the solution via an L^—contraction property; notice that this contraction is valid for a general 
open set 0 with boundary. 

Let us write the weak formulation of (3.2): for every smooth test function cf : [0,T] x fl —>• M 
with V(?i> • h = 0 on [0, T] x dLl we have 


(3 


. 3 ) r [ [pd, 

Jo Jn 


k4> + {p + p)A(f + pu ■ Vcp] dx dt + / pQ(x)(j){0,x) dx = / p{T,x)(j){T,x) dx. 


In 


By density arguments the above formulation holds for f G T]; L^(f2))nL^([0, T]]H‘^{Ll)). 

Now, let us consider two solutions {p^,p^) and {p^,p^) of Problem 3.2, with initial conditions 
respectively pj and Pq. Writing the weak formulation (3.3) for both of them and taking the 
difference we obtain 

Z{(l), T) = 0)+ [ [ [(p^ - p‘^)dt(f) + (p^ - p^ + p^ - p^) A(/> + (p^ - p2)u • Vcf] dx dt, 

Jo Jn 

where T{(j),t) = / (j){t, x)[p^{t, x) — p‘^{t,x)] dx. We introduce the following quantities 
Jn 


A := 


p^-p^ 


and B := 


p^ — p^ 


(pi - p2) + (pi - p2) ■ (pi - p2) (^pl - p2y 

Note that 0 < ^ < 1 and 0 < B < 1 a.e. in [0,T] x and A + B = 1. To be consistent with 
these bounds, we set ^ = 0 when pi = p^, even if pi = p^ and B = 0 when pi = p^, even if 
p^ = p^. With these notations the weak formulation for the difference gives 

(3.5) X((/>, r) = X((/), 0) + f f {{p^ - p'^) + {p^ - p‘^))[Adt4> +{A +B)A4> +Au-\/(f] dxdt. 

Jo Jn 

For a smooth function 0 : ^ M such that \0\ < 1, let us consider the dual problem 

Adt4> + (^ + B)A(j) + Au • V(p = 0, in [0, T) x Ll, 


(3.6) 


V0 • h = 0 on [0, T] X dLl, 4>{T, ■) = 0 a.e. in Ll. 


Let us remark that if we are able to find a (reasonably regular) solution (j) for this problem 
for any 6 smooth then, using the maximum principle and then optimizing in 9, we would get 
an L1—contraction result for p: in particular we get uniqueness for the initial value problem for 
p and hence also for p (as done in the end of the proof of Theorem 2.4). However, since the 
coefficients in (3.6) are not regular, we study a regularized problem. For e > 0 let us consider 
A/r, B^, Ue to be smooth approximations of A, B and u such that 

(3.7) M — ^£||L>'([o,r]xr2 ) + 11^ “-S£||l’'([ 0 ,T]xC) + l|u — U£||j;^rQo,T]xO) < ^) ) e<^e,Be<l, 

for some 1 < r < +oo, the value of which to be chosen later. The regularized problem reads as 
follows 


(3.8) 


dtcpe + (1 + BefAyAcj)^ + Ue • V(/>£ = 0, in (0, T) x 12, 

V(^£ • h = 0 a.e. on [0, T] x dLl, 4>e{T, ■) = 9 in Ll. 
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For all e > 0 the above problem is uniformly parabolic and B^jA^ is continuous and positive. 
Moreover 0 is smooth, thus by classical results (see for instance [20, 19]) the problem has a 
(unique) solution (jie G <^^([0, F] x O). In particular can be used as test function in (3.3). In 
the followings we shall use some standard uniform estimates (in e) on given in Lemma 3.2. 
In particular, using as test function in (3.5) one has 

fT r 

x(0„r)-x((/.„o) = 


'0 Jn 

/7' 

Jo Jn 

rT 


{p^ - +p^ - 7) [Adtpe + (^ + B)Aps + • S/pe] dx dt 

{p^ - p^ +p^ - 7) [Adtpe + (^ + B)Ap^ + • Vps] dx dt 

[ [ {p^ - +P^ - 7)^ [dtps + (1 + Be/A^)Ape + Ue ■ Vpe] dx dt 

Jo Jo. 

= [ [ {p^ - p‘^+P^-p‘^)iBs/Ae)iAe-A)Apedxdt 

Jo Jn 

+ f I {p^ - p"^ + p^ - p‘^){B - Be)Apedxdt 

Jo Jn 

+ [ [ {p^ - p"^+p^-p‘^)A{ue-u)-Vpedxdt 

Jo Jn 


>0 Jn 

■- /I _g/2 j3 


Let us show that |/j| —>■ 0, |/|| ^0 and |/|| ^ 0 as e —>■ 0. More precisely, first let us recall 
that 0 < < 1 a.e. in [0,T] x 12, hence p^iP^ G L°°([0,T] x ft). On the other hand 

p^,p^ G L^([0, T]; Lt^(f2)) and by Corollary 1.2 we have that 

[ |Vpipdx< [ jutpdt, 

Jn Jn 

for almost every t G [0, T]. This implies that (since u is bounded) 

ess - suptg[o,r]l|Vpi||L2(f^) < C. 

In addition, p*’s being pressures one has \{pl = 0}| > \{pl < 1}| > |n| — 1 > 0 for a.e. t G [0, T], 
and so by a suitable version of Poincare’s inequality (since pi vanishes on a set of positive 
Lebesgue measure) one obtains p* G L°°([0, T]; In particular, by the Sobolev embedding 

this implies that p* G L°°([0, Tj; L'^(f2)), i G {1, 2} for all 1 < < 2*. Let us fix g G (2, 2*), where 

2* = 2d/{d - 2) if d > 3 and 2* = oo if d = 2. 

This implies the following estimates 

l7l < ||7 “ 7 ||l°°([0,T]xC) • WiBe/AeY^'^iJ^ - Ae)\\L^(l0,T]xn) ' II (-SeMe)^'^^ A(/)£ ||^2([o,r] xO) 

+ r 11 / -7||l.(0) • \\{Be/Ae)^/\A - Ae)hr^n) • \\{Be/Ae)^/^APeh^n) dt 
Jo 

< C{l/ef/^e 

+ 11/ -/||l-(L.) • \\{Be/Aefl\A- Ae)\\L2^Lr) ' \\{Be/Aef^^APeU^i^^p 

< ^0, as e ^ 0 


and similarly 

|/| < 11/ - P^||l,°°([0,T]xO) • \\{Ae/Be)^^'^{B - i^e) ||l2([o,t] xO) ' \\{Be / ^Pe\\L'^{[0,T\xn) 

+ 11 / -/||l^(L.) ■ \\{Ae/Bef/\B - Be)h^^Lr) ' || (B.M.)'/'||^ 2 (^ 2 ) 

<C{l/e)^/^e = Ce^/^ ^0, as e —>■ 0; 
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finally we have 

I < \\p^ — • ||u — Ue||2,i(i:,r) • ||V(/)£||x, 00 (^ 2 ) < Ce ^ 0, 

where r > 1 is the exponent such that —I-1— = 1, i.e. r = — ^ (recall that q G (2,2*). We 

2 r q O' — 2 

used the approximations (3.7) with this specific r. Hence we obtained that 

= 2 :( 0 £,O) + 0 ( 1 ); 

now we use the fact that we fixed the final condition 9 for all e, and also that |(/>e| < 1, that is 
implied by the strong maximum principle (see for example Theorem 2.9 in [24]); in particular 
we can write 

[ 9{x){p^{T,x) - p‘^{T,x)) dx = [ (I)si0,x){pl - pl)dx + o{l) < [ \pl - pl\dx + o{l). 

Jq Jq Jq, 


Letting e ^ 0 we finally find / 0{prp — pjPj dx < ||Po “ PoIIrI) optimizing in |0| < 1, we 

Jn 

proved the contraction. In particular this implies that p^ = p^ a.e. in [0,T] x H whenever 
Po ~ Po- Then one obtains also that = p^ a.e. in [0,T] x fl, as in the end of the proof of 
Theorem 2.4. The result follows. 

□ 


Lemma 3.2. Let he a solution of (3.8). Then there exists a constant 

C = C{T,\\u\\L^,\\V9h2^^))>0 
such that we have the following estimates, uniformly in e > 0 .' 

(i) sup \\V(peit)\\L 2 (n) <C-, 
te[o,T] 

(ii) \\{Bs/Ae)^A4)s\\L^([o^T]xn) AC] 

(hi) ||A0e||£,2Qo,T]xO) < C. 

Proof. Let us multiply the hrst equation in (3.8) by A(j)^ and integrate over [t,T] x H for 0 < 
t < T. We obtain 

2 ll^‘?^e(^)llL 2 (C 2 ) + [ f (^ + Be/A,.)\A4>s\‘^ dxdt 

(3.9) [ [V^pdx- [ f Ue -VfieAcpsdxdt 

2 Jn Jt Jn 

Hence by Young’s inequality we have 

Iw^Mmh^n) \\^Ms)\\h^n) d^ + ^l|V0||i2(^) 

<C'+^^^||V</-,(s)||i,(^)ds 

where the term in |A(/>£p has been absorbed by the left hand side, and 0 < 5 < 2/||u||ioo is a 
fixed constant and the constant C > 0 is depending just on ||V0 ||j;^2(q) and ||u||ioo. Hence by 
Gronwall’s inequality we obtain 

which implies in particular that sup ||V(/>£(t)||i 2 (Q) < C. Thus (i) follows. 

tG[0,T] 
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On the other hand choosing 6 := 2/||u||ioo in Young’s inequality used in (3.9) and using (i), 
we obtain 

f f {B,/Ae)\A(P,\^dxdt<C 
J t J Q, 

hence ||(Se/A£)^A(/)£||i2([o,T]xO) < C, and thus (ii) follows. 

By (3.9), (i) and (ii) easily imply (hi). □ 

4. About the L^—contraction in the first order case 

In the previous section we proved uniqueness in the second order case using an contraction 
result. We expect this result to be true also in the first order case; however we expect the treat¬ 
ing of this L^—contraction problem in the most general framework (let us say when u £ W^’^) 
to be difficult, since this approach should “include” (in the proof) also the well-posedness in the 
Di Perna-Lions theory as a special case. In fact, when one has V • u = 0 and if one starts with 
Pq then this condition is preserved without adding the pressure term. 

However in the case treated in Section 2 (i.e. with a monotone velocity field) we can try 
to sketch a proof of the contraction result using the uniqueness already proved: let us 
approximate the solution discretizing in time using the splitting discrete scheme “continuity 
equation -|- Wasserstein projection onto the set {p < 1}” (similarly to the scheme introduced in 
[31]). Since the vector field is monotone there exists a unique flow, that implies also that there 
exists a unique solution to the initial value problem dtpt + V • {yitPt) = 0, pt^ = p; in particular 
one can define a function 'I'^J(/i) := pt^, that will satisfy also the semi-group rule 'kf o'k* = 'I'| . 
Let T be a time step and define recursively Pq = po and then 

T _ [ if ^ is even 

\V[pn) if n is odd. 

Then the distance is preserved through the continuity equation step while it decreases after 
the projection thanks to Lemma 4.1. So in particular the contraction is true in the discrete 
scheme; once we have that the scheme converges as r —)• 0 to our equation, the uniqueness 
result gives that this property is preserved in the limit. In order to guarantee convergence 
(see [10] for general conditions) the crucial quantity to estimate is IL 2 (Tj“'"'”(p), p); using the 
Benamou-Brenier formula, two conditions that guarantee a good estimate are: 

• u G L^([0, T]; L°°(M'’*)): in this case we would have W|('I'j“'"'”(p), p) < r Uu^jj^oo ds; 

• (V • u)_ G L^([0, T]; L°°(M'^)): in this case we would have [['Ll(p)jj 2 ,oo < Cjjpjjioo for 

some universal C and in particular we have {p)-, p) < CII/jIIloot ||usll ^2 ds. 

We believe that this general scheme (uniqueness in the Wasserstein framework and approx¬ 
imation with contractive time discrete approximations) could be adapted also with some 
other convection terms 4>. 

Lemma 4.1. Let us consider the projection operator V : —)■ 

(4.1) 'P{p) = argmin{lL|(p, r/) : p G p < 1}. 

Then, when pi,p 2 are probability densities we have \\V{pi) — 'P{p 2 )\\l^ < [[pi — P21 Ili- 

Proof. First we see that formula (4.1) let us extend V to measures in ./#+(M'’*) with mass less 
than jHj. In this context a monotonicity property is true (see for example Theorem 5.1 in [1]), 
that is we have V{p) < 'P(p) almost everywhere if p < p. Now we can derive the contraction: 
let us denote p = min{pi,p 2 }. First of all we have that 

|pi - P 2 I = max{pi,p2} - min{pi,p2} = (pi -b P2 - min{pi,p2}) - min{pi,p2}; 
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in particular ||/9i — P 2 ||li = 2 — 2 f pdx. By the monotonicity we have 'P{pi) > V{p) and in 
particular we have m.m{V{pi){p 2 )} > ^{p) and so 

\\'P{Pi)-'P{P 2 )\\l^ = 2-2 J mm{V{pi),V{p 2 )}dx < 2-2 J V{p)dx = 2-2 J pdx = ||pi-/32||lB 
which proves the claim. □ 
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